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Abstract 

The role of nuclear spins in decoherence and dephasing of a solid state 

phase qubit is investigated. Both effects of static spin environment and spin 

polarization fluctuations in time are considered on the basis of non-Markovian 

Langevin-Bloch equations. We find conditions when coupling of a phase qubit 

to a bath of nuclear spins does not impair coherent evolution of the qubit. 
76.20.+q, 75.45.+j, 03.67.Lx 
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I. INTRODUCTION 

An analysis of decoherence and dephasing in quantum systems interacting with an en- 
vironment has attracted a great attention in the recent years. This problem is of particular 
importance for a physical implementation of quantum computers, and, especially, for solid- 
state designs of qubits Electronic or nuclear spins embedded into solid state materials 
(see papers by D.G. Cory et.ai, G. Burkard et.ai, A. Imamoglu, and by B.E. Kane in Ref. 

provide a natural way for a practical realization of quantum bits. However, difficul- 
ties associated with a manipulation of individual spins render the spin-based design quite 
challenging. Superconducting qubits are more controllable, besides that they can be 
produced by means of well-developed litographic methods and offer advantages of scalability 
and flexibility. Unfortunately, operating states of a superconducting charge qubit extensively 
studied in the last time are highly coupled to environment degrees of freedom that leads 
to the high level of decoherence in the system. To get around this obstacle several models of 
phase superconducting qubits were proposed 0-0. These qubits utilizing superconducting 
phase as a dynamical variable with two macroscopically distinguishable values are effectively 
decoupled from charge degrees of freedom. In the general case the characteristic property of 
the phase qubit is the presence of a small magnetic flux which has two opposite directions 
in the different states of the qubit. An interaction of the magnetic field created by the 
fiux with a bath of nuclear spins will be an important mechanism resulting in destroying 
a quantum superposition of qubit's states and in statistical dispersion of characteristics of 
different qubits. In this work we develop a mathematical model describing the contribution 
of nuclear spins to decoherence and inhomogeneous broadening of phase qubits and deter- 
mine conditions when nuclear spins have only a slight effect on the operation of the qubit. 
The problem of the spin-bath mediated decoherence in phase qubits has been previously 
discussed in Refs. P-p^. In papers particular attention has been given to the case of 
static spin environment when the nuclear or paramagnetic spins have been subjected to the 
action of many magnetic fields, including the magnetic field from the qubit together with an 
external magnetic field and dipolar fields from other spins. By means of instanton technique 
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Prokof 'ev and Stamp have obtained, in particular, some estimations for the dephasing 
rates of the phase qubit in the static spin environment. It should be noted that the simple 
model considered in our paper allows us not only to get the same estimations but also to 
describe in detail a time evolution of the qubit 's variables, and, especially, to determine the 
spectrum of all possible frequencies of quantum beats in the ensemble of qubits. This is a 
must for a correct analysis of logical operations in the ensemble of qubits. We suppose here 
that the magnetic field created by the qubit dominates over other magnetic fields acting 
on the nuclear spins. This assumption would hold for the qubit 's design based on small 
grain-boundary junctions between d-wave superconductors . 

Besides the static case, in the present paper we consider a contribution of spin-lattice 
relaxation to the decoherence rate of the phase qubit. This phenomenon which is completely 
ignored in the above-mentioned paper leads to true decoherence as distinct from 



dephasing caused by the static spins. A qualitative estimation of effects of time-dependent 
spin fluctuations on the superconducting persistent current qubit has been performed in Ref. 
pi| . In the paper |TD| nuclear spins are supposed to be in the state of thermal equilibrium 
despite the fact that the qubit 's magnetic field oscillates with a high frequency determined by 
splitting of two lowest energy states. Under such conditions an application of the fluctuation- 
dissipation theorem for the calculation of the qubit 's decoherence rate as it has been done 
in Ref. 0] casts some doubt. A longitudinal relaxation of nuclear spins takes long enough 
time Ti, however, many nuclear spins will contribute to qubit's decoherence, so that the total 
decoherence rate of the qubit will be proportional not only to the damping rate of a single 
nuclear spin but also to the number of spins in the qubit's area. It follows herefrom that 
the mechanism of qubit's decoherence due to spin-lattice relaxation is competitive with other 
mechanisms and, at least, can not be rejected without more circumstantial consideration. 

The present paper is organized as follows. The next section is devoted to the formulation 
of the problem. In Sec. Ill we analyze inhomogeneous broadening in the ensemble of qubits 
caused by a static spin environment. Decoherence rates of the qubit interacting with a spin 
polarization fluctuating in time are calculated in Sec. IV. 



II. MODEL 

As mentioned above, there is a spontaneous flux ($ ~ 10^^ — 10^^ $o) concentrated 
in the central part of the phase qubit and having opposite directions in two degenerate 
equilibrium states iyjo |@H31- This flux creates the corresponding magnetic field Bo(r) which 
also has two directions: along and opposite of z-axis, ±i?o(r) at different values of the qubit's 
phase. The existence of two-well potential for the phase difference in the asymmetric grain 
boundary Josephson junctions in d-wave cuprates directly follows from the experimental 
results obtained in Ref. [0. A doubly-degenerate ground state of such a junction is split into 



the symmetric and antisymmetric states with A being the energy difference between them. 
The qubit's two-level system is conveniently described by the Pauli's matrices Tj.,Ty,Tz, so 
that the tunneling Hamiltonian is given by the expression: Hq = {A/2)tz. The eigenvalues 
±1 of the matrix corresponds to the positioning of the qubit in the right well {(fo) or 
in the left potential well (— y^o), whereas the eigenvalues of the matrix corresponds the 
energy levels ±(A/2) of the antibonding (antisymmetric) and bonding (symmetric) states. 
The magnetic field produced by the qubit can be also represented as a matrix: Bq = TxBo{r) 
with JB{r)dS = ^. 

Nuclear spins in the area of the qubit are subjected to the action of this magnetic field, 
and their energy operator looks as follows ||T0[| 



Hs = -T^^g^iNiB{ri){ai)z (1) 

i 

where is the nuclear magneton, (crj)^ is z-projection of the Pauli's matrix describing 
i — th nuclear spin, B{ri) is the magnetic field induced by the qubit in the point where 
i — th nuclear spin is located. Here we take a sum over all nuclear spins placed in the qubit's 
area. The effect of any external magnetic field on the nuclear spins is supposed to be small 
compared to the action of the qubit's magnetic field. It might be well to point out that the 
Hamiltonian (1) can describe an interaction of the qubit not only with nuclear spins but 
with any localized magnetic moments as well, and, in particular, with an electronic spin 1/2 



degree of freedom located in the Cu02 planes of the YBCO superconductor [|T4|. In this 



case we have to extend the sum in Eq.(l) to the electronic spins with the substitution of 
electronic magnetons ^3% for ^^Ni1 and the consideration cr, as electronic Pauli's matrices. 

To study an internal dynamics and fluctuations in the system of nuclear spins we sup- 
pose that i — th nuclear spin is coupled to the dissipative environment with variables 
{Qx{f)^Qy{t),Qz{t)}i. This coupling has the form 

HN = -Y.^^■^^■ (2) 



For example, a heat bath of acoustic phonons is described by the variable [|I5 



Here h^^h^^ are creation- annihilation operators of acoustic phonons with wavevectors q and 
the sound velocity c^; G is a constant of spin-phonon coupling, p and V are the density and 
the volume of the crystal, respectively. The spin-phonon interaction (2), while unavoidable, 
is weak enough, especially for the nuclear spins having a small magnetic moment. As a 
result, the directions of nuclear spins do not conserve exactly. But due to the weakness 
of spin-phonon coupling we can suppose in the first approximation that nuclear spins have 
fixed directions which, however, vary from one qubit to another. In the part III we investi- 
gate inhomogeneous broadening or dephasing in the ensemble of qubits associated with this 
variation. Time-dependent spin fluctuations due to the spin-phonon interaction lead to the 
true decoherence of the qubit. The section IV is devoted to the study of this phenomenon. 
It is of interest that the spin-lattice relaxation for the copper, oxygen and yttrium nuclei is 
usually considered as a result of the interaction between nuclear spin and electronic spins 
localized in the CUO2 planes [1^. Besides that, the electronic spins contribute directly to 



the decoherence of the qubit with the polarization, P, containing the time-dependent part 
only. 

The total Hamiltonian of our system will look like: 

H=^r,- Pr. - ■ + Hb. (3) 



Here 

P-T.^{^i)^ (4) 

is the operator of Zeeman energy of nuclear spins in the fixed qubit's magnetic field, Aj = 
2giif^iB{ri) is the Zeeman frequency of i— th nucleus, and Hb is the Hamiltonian of the free 
heat bath. In the case of acoustic phonons the Hamiltonian Hb will look like 

^B = E^c,|q|(6+6q + l/2) 

We suppose that the free heat bath (without the interaction with the spins) is in thermal 
equilibrium with a temperature T. The operator P (4) is proportional to the spin polar- 
ization in z-direction. Because of this, we will mention this operator later just as a "spin 
polarization" . 

It follows from the Hamiltonian (3) that the time evolution of the qubit's quasispin 
described by the matrices Tx,Ty, is governed by the equations 

fy = + 2r,P, 

= -2ryP. (5) 

Direct and indirect interactions between different nuclear spins are neglected in this 
model. In view of a great number of nuclear spins in the qubit's area an influence of these 
spins on the qubit's dephasing can be quite considerable. Fortunately, the main part of the 
spin polarization, Pq, is almost static and determined by the initial spin configuration of 
the given qubit. This part differs in different qubits, as well as in the same qubit but taken 
at the time of the next logical operation. These fluctuations result in statistical dispersion 
of properties of qubits followed by the impact on the entanglement of qubits. In the next 
section we will analyze a contribution of the static spin polarization without resorting to the 
perturbation theory. Besides the static part the spin polarization will include a small part 
fluctuating in time. This part P{t) — P{t) — Pq is supposed to be here due to the spin-lattice 
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interaction and results in destroying of quantum coherence of each individual qubit. 



III. INHOMOGENEOUS BROADENING 

In this section we consider effects of the static spin environment {P — Pq) on the time 
evolution of the qubit's ensemble without taking into account the component fluctuating in 
time. The interaction between the qubit degrees of freedom and the static nuclear spins is 
not supposed to be weak. Each of the qubits has a specific distribution of nuclear spins. 
We have to average qubit's operators over all configurations of nuclear spins to obtain an 
information about a time evolution of the qubit's ensemble, for example, an information 
about frequencies of quantum beats (tunneling oscillations) of different qubits. Our goal is 
to find conditions when the frequency dispersion of qubit's ensemble due to the interaction 
with nuclear spins will not influence on a performance of logical operations, especially, on 
the entanglement of different qubits. 

In this case the time evolution of the qubit's operators is described by the equations (5) 
with the time-independent spin polarization P — Pq. A direct solution of the system (5) of 
differential equations (with P — Pq) gives us the formula 



1- ^(1-cosQi) 



T.(0) 



A^^T,(0) - ^(1 - cosni)T,(0) (6) 

that represents oscillations of qubit variables with a frequency = + 4Pq in a specific 
spin environment characterized by the parameter Pq- However, this parameter is changed 
when we proceed to the next qubit. To study a variety of possible evolutions we have to 
average the expression (6) over an ensemble of all spin configurations. It has been done in 
the Appendix A. As a result, the time evolution of the qubit's variable averaged over nuclear 
spin configurations is given by the following expression 



ra=N 
m=-N 



^ cos(t^4Agm2 + A2) 



4A2m2 + A2 



(r.(0)) 



^ „ sinitJA/^lw? + A2) 
E C^N-^ 'l'^ \ r,m. (7) 

m=-N JAAlm^ + A2 



m=N 
1-2N \ " r^N-m 



In the case of zero Zeeman frequency (Aq = 0) the phase differences of all identical 
qubits in our ensemble will oscillate with the same frequencies equal to tunnel splitting A. 
Otherwise, the frequencies of quantum beats of different qubits can take the spectrum of 
values ujjn — \J 4 Apm^ + A^ with m — 0, 1, .., A^, because each of these qubits has the own 
quasistatic distribution of nuclear spins (up or down), and, therefore, the own energy of 
interaction with these spins. Changing an orientation of a single nuclear spin, e.g., from up 
to down, will change the interaction energy between this spin and the qubit by the value of 
Aq. This fact results in the discrete character of qubit's frequencies. During the time of one 
logical operation the frequency of the given qubit is equal to one of the frequencies But in 
the time of performing of the next operation separated from the first operation by the interval 
that is in excess of the spin-lattice relaxation time Ti this qubit can be characterized by 
another frequency from the set {o;^}. The adjacent qubits can also have different frequencies 
that will affect entanglement processes. The dispersion of qubit's freqiencies is described by 
the binomial coefficients 



{N-m)\{N + m)\ 



times 2 . At large N and m <^ N this function is approximated by the Gaussian distri- 
bution: 

2-2^C2^-'" ~ exp(-mVAr) 



so that the probable values of m are in the range from zero to v A^. The maximum frequency 
of quantum beats is of order ^4AqA^ -|- A^. In the practically interesting case of weak 
coupling between nuclear spins and the qubit we obtain for the inhomogeneous broadening 
of the qubit's Une 

{AujUak = 2N^, (8) 
whereas for the strong coupling we have 
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In the framework of instanon technique Prokof'ev and Stamp |^ have considered more 
general model of qubit's dephasing due to static nuclear spins, but at the final stage they 
have obtained the same estimations as ones in Eq.(8),(9). This points to the ability of 
our simple mathematical model to describe properly the above-mentioned mechanism of 
qubit's dephasing. It should be emphasized that besides of estimations our model gives the 
mathematical expression for the time evolution of the qubit's variable {Tx{t)) (7) that is of 
great importance for an analysis of logical operations in the ensemble of qubits. 

The static nuclear spin environment does not affect the time evolution of the qubit 
and, therefore, the performance of logical operations, if the number of nuclear spins, 2N, 
interacting with the qubit is less than the ratio of tunnel splitting to Zeemann splitting 
squared: 



2N<^, (10) 



For the sake of estimation we consider the qubit's model with the linear size of the central 
area / ~ 10~^cm. This size seems to be too small for designs considered in 0,0, but it is 
very reasonable for qubits based on Josephson junctions between d-wave superconductors 
Pj6| JTl|JT^ . As shown in [|1^] for this design the magnetic field produced by circulating spon- 



taneous currents can be localized in the area with one size of about lO^o (^o is the coherence 
length, ^0 ~ lO^^cm) and another size of about the width of the junction (~ 0.2yum). Such 
small junctions are necessary also to eliminate experimentally defects and other irregulari- 
ties. Even with both sizes of order 0.2/xm and with the magnetic flux created by the qubit 
of order $ ~ 10^^<I>0) '^'o = hc/2e we obtain that the qubit's magnetic field can reach the 
value 10 G comparable to the dipole fields between the nuclear spins. Different components 
of Y Ba2Cu'iO'j have magnetic moments in the range from — 0.238(®^y) to 3.075(®^C'u) of 
nuclear magneton. Because of this, for the Zeeman frequency of one nucleus in the qubit's 
magnetic field we get the estimation Aj ^ 10'^ — 10^(l/s), so Aj/A ^ 10^^ — 10^"^ ^ 1, if 
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A ^ 10^(l/s). As follows from Eq.(lO), the number of nuclear spins coupled to the qubit 
should be less than 10^ — 10^ if we want to get rid of effects of above-mentioned inhomoge- 
neous broadening on the qubit's evolution. The magnetic field of the qubit can not penetrate 
into the superconductors at the depth in excess of the penetration length A. Only nuclear 
spins located inside of the volume with third dimension smaller than the penetration length 
will be subjected to the action of the qubit's magnetic field. For a Y Ba2Cu^0j sample 
with cell sizes of order (3,3,11)A there is almost 2N ~ 10^ nuclear spins in the qubit if 
the penetration length A ~ 10~^cm. This number of nuclear spins is close enough to the 
above-mentioned requirements. With small decreasing the number of spins (or decreasing 
the ratio Aq/A^) the statistical dispersion of frequencies of different qubits will be of little 
importance and the operation of the qubit under discussion will not be distorted by this 
mechanism of inhomogeneous broadening. It should be mentioned that the qubit's magnetic 
field and, therefore, the Zeeman frequency of a single nucleus decreases with increasing the 
size / of the qubit: Aq ~ *&/^^, whereas the number of nuclear spins, A^, is proportional to 
the qubit's area, N P x \, provided the qubit's magnetic flux, $, and the penetration 
length. A, remain constant. Then, the inhomogeneous broadening of the line (8) is inversely 
proportional to the area of the qubit: {Au)weak ~ 1/^^, and we can expect that for larger size 
qubits studied in Refs. 0,|lOl the effects of nuclear spins are of lesser importance. However, 
in this case the magnetic field created by the qubit can be less than the dipole field between 
the nuclear spins, and we need the generalization of the formulae (8), (9), (10) to describe 
the situation correctly. 



IV. DECOHERENCE OF PHASE QUBIT IN THE SPIN BATH 

Now we analyze an effect of a nuclear spin polarization fluctuating in time on decoher- 
ence processes in the superconducting qubit. The contribution of electronic spins localized 
in Cu02 planes of the YBCO superconductor to the qubit's decoherence rate can be 



calculated by the same approach. We assume here that the concentration of nuclear spins are 
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sufficiently small, so that the condition (10) is fulffiled and we can neglect an effect of static 
spin polarization on the evolution of the qubit: Pq = 0, so that fluctuations of the spin po- 
larizations are completely due to coupling of nuclear spins to the heat bath, e.g. to phonons. 
Now a contribution of the fluctuating component to the total Hamiltonian (3) of the system 
is supposed to be small compared to the tunneling term {A/2)tz. Here we consider also 
weak coupling of nuclear spins to a dissipative environment, for example, to a phonon heat 
bath, as a source of fluctuations of nuclear spin polarization. This coupling results in the 
energy dissipation and thermalization of the spin system. In the weak-damping approxima- 
tion when the energy of the interaction between spins and the heat bath is much less than 
the Zeeman frequency A, we can find a dissipative evolution of nuclear spin operators a^j in 
the ffipping magnetic field TxBziri) created by the qubit and calculate a correlation function 



Here brackets (..) mean averaging over the state of the equilibrium heat bath. It should be 
noted that the equations (5) descriptive of the qubit 's evolution incorporate the total spin 
polarization P = P + (P). However, the mean part of the polarization (P) leads only to 
small corrections for the tunneling frequency A, and, therefore, it can be dropped out. At 
the same time the fluctuating part of the spin polarization P causes the destruction of the 
qubit 's coherence for the time interval which is equal to the relaxation time of the averaged 
X— projection (r^.) of the qubit's quasispin. The evolution of matrices t^, Ty, is governed 
by the equations (5) where P = P{t). The fluctuating part of spin operator P{t) includes 
contributions from many different spins and can be described approximately by Gaussian 
statistics. With this property in mind the mean value of the product of the qubit's matrix 
and the spin polarization, say, (r^P), can be expressed as |]T^|-|T^] 



{[P{t),P{t')]^) = {P{t)P{t')±P{t')P{t)) 



with the polarization P{t) being determined by Eq.(4): 
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Here 

Mp{t,t,) = {^[P{t),P{h)U) (11) 

is the correlation function of the spin bath in the presence of flipping magnetic field created 
by the qubit, [A, B]± = AB ± BA,Mp{t,ti) = Mp{t,ti)e{t - ti) with e{t - ti) being the 
Heaviside step function. Averaging Eq.(5) over fluctuations of the Gaussian spin variable 
P{t), i.e. eventually over the heat bath fluctuations, gives the non-Markovian equations 
|16H18| for the mean qubit 's variables 

{fy) = A(r,) + 2 j dhMpit,ti){t[T,{t),r^{h)]^), 

Tz = -2 j dtiMp{tM)^[ry{t),T^{ti)].). (12) 

As mentioned above time-dependent fluctuations of spin polarizations are supposed to be 
weak. Because of this, to calculate the commutators in Eqs.(12) we can use the equations 

+ AV, = 

for the free oscillations of qubit's variables without any spin polarization. Then, the com- 
mutators are expressed in terms of qubit's matrices, for example, 

i[Tz{t), r^.(ti)]= - 2Ty{t) cos A(t - ti) + 2r^(t) sin A(t -ti). 

Ignoring the frequency shift of the qubit due to spin fluctuations we obtain the simple 
equation 

(f,.) + 2r(f,.) + A2(r,.) = (13) 
describing the relaxation of qubit's oscillations between two wells with the relaxation rate 

r = 5p(A), (14) 

where 
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/oo 
dre^'^^Mpir) (15) 
-oo 

is the spectral density of spin fluctuations. The relaxation of the qubit matrix Tx to the 
equilibrium value is described by the expression : 

(r.(t)) = [(r,(0)) cos(At) - (r,(0)) sin(At)] 6"^*. (16) 

The coefficient F (14) represents also the rate of qubit 's decoherence in the spin envi- 
ronment. To find this decoherence rate we have to calculate the spectral function Sp{u!) 
(15) and the correlation function Mp(t,ti) (11). Fluctuations of different spins are sup- 
posed to be independent. In view of this fact together with Eq.(4) the correlator Mp{t,ti) 
is proportional to the sum of correlation functions of ^^-projections of nuclear spins: 

Note that only 2;— projection of the nuclear spins is coupled to the qubit because the qubit's 
magnetic field is always directed parallel to z-axis while its magnitude oscillates with time. 
The calculation of the correlation function of i-th spin ((l/2)[((Tj)^(t), ((Ti)^(ti)]+) in the 
presence of hipping magnetic field and a heat bath together with finding the spectral density 
of spin fiuctuations Sp{u!) (15) has been performed in the Appendix B. 

With the formulas (14),(B21) we obtain the following expression for the decoherence rate 
(or the inverse decoherence time t^^ ) of the phase qubit 

r _ 1 - 1 VA2 7»(A) 

^-^-2V''^A^ + 7f(A)- ^''^ 

To obtain the expression (18) we have supposed that the decoherence rate, F, is much 
less than the qubit's tunneling frequency, A, and that the the linewidth of i-th nuclear spin, 
7j, is less than the corresponding frequency A^: F <^ A,7j -C Aj. 

In view of the fact that the relaxation rate of i-th nuclear spin, 7i(A), or, in other words, 
an inverse longitudinal relaxation time (1/Ti)j of i-th spin, is much less than the qubit's 
frequency A : (1/Ti)j = 7i(A) <^ A, we find that the qubit's decoherence rate is proportional 
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to the sum of spin line widths 7^ due to spin couphng to a heat bath taken with the weights 
(A^/A^) : 

This formula takes into consideration oscillations of the qubit's magnetic field with the tun- 
neling frequency A. Usually Zeeman splitting Aj of i-th nuclear spin in the qubit's magnetic 
field as well as in the magnetic field created by other nuclear spins is much less than the 
frequency of quantum beats A : Aj ^ A. Because of this, the term with m = gives the 
major contribution to the sum in Eq.(19), and in this approximation the qubit's decoherence 
rate will look like: 

l = 2,"(A)coth(A)i:§. (20) 

Here x"{^) is the imaginary part of the susceptibility of the heat bath with the temperature 
T. The interaction of nuclear spins with this heat bath, for example, with lattice vibrations, 
provides a mechanism of dissipation and thermalization in the spin system. If the heat 
bath temperature is in excess of the tunnel splitting: T ^ A, the qubit's decoherence rate 
increases linearly with temperature: 

The nuclear spin-lattice relaxation time is usually measured at the field of order lO^Gs 
when the Zeeman frequency for one of our nucleus is 10^(l/s). At this field and at the 
temperature T ~ IK the spin-lattice relaxation rate is of order or less than 10(l/s) ||19|| . 
According to results by Morr and Wortis pO[ the spin-lattice relaxation rate is approximately 



proportional to the applied magnetic field squared. Our decoherence time, (19)-(21), is 
determined by the spin-lattice relaxation rate 7i(A) taken at the frequency A. For the 
qubit's model mentioned in the end of Sec. Ill we have A ^ 10^(l/s), that is more than the 
usual frequency 10^(l/s) of relaxation time measurements by a factor of 10. So, the spin- 
lattice relaxation rate corresponding to our quantum beats frequency A will be estimated as 
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(1/71)^ = -fi(A) ~ 10^(1/5). For the same ySaaCiigOy sample as in Sec.III with 2N ~ 10^ 
nuclear spins the qubit's decoherence rate due to magnetic coupling to the nuclear spins is 

- ~ 10^(l/s). 

This decoherence time is much more than the period of quantum beats oscillations: r^A ~ 

10^ 3> 1, so that the necessary condition 

^tunneling ^ '^gate ^ '^decoherence 

or 

holds if the nuclear spins are considered as the main decoherence mechanism. Nevertheless, 
as is evident from the foregoing, the mechanism of the qubit's decoherence related to the 
spin-lattice relaxation can not be considered as improbable, especially as this mechanism 
causes true decoherence of a single qubit in distinction to dephasing occuring in the en- 
semble of qubits. It is necessary to emphasize that the role of localized magnetic moments 
related to the electronic shells in the superconductor has to be investigated as well for the 
correct estimation of the decoherence rate r^^. To estimate the contribution of electronic 
spins belonging to the Cu and O ions to the decoherence of the qubit we can assume that 
the ratio between the electronic spin magnetic moment and the nuclear magneton is about 
10^, so that for parameters mentioned at the end of the part III, the Zeeman sphtting of the 
electron in the qubit's magnetic field is of order Aj = Aq ~ 10^(1/ s). The spin line width, 
7j (B19), increases more significantly because it is proportional to the electronic magneton 
squared: 7, ~ 10^ — 10^(l/s), and reaches the values of Aj and A. It is evident that for the 
correct analysis of this limit we have to generalize the formula (18), however, this task is 
beyond the scope of the our paper devoted to the effects of nuclear spins. Electronic spin 
magnetic moments contribute, in particular, to the spin-lattice relaxation of nuclear spins, 
and, by doing so, indirectly affects the decoherence rate of the qubit. It is of interest that 
the decoherence rate (18) grows with the value of the spin magnetic moment showing the 
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tendency for a saturation. Because of this, we can expect that the locahzed electronic spins 
will not destroy the working conditions for the phase qubit. 



CONCLUSIONS. 

We have considered a contribution of nuclear spins to the dephasing and decoherence of 
a phase qubit. An effect of the static spin environment has been analyzed without resorting 
to the perturbation theory. For this case we have found a condition (10) when dephasing 
or inhomogeneous broadening in an ensemble of qubits has no effect on the operation of a 
quantum computer. Decoherence rates (19), (20), (21) of the qubit weakly coupled to the heat 
bath of nuclear spins have been calculated as well; in so doing we have taken into account 
fluctuations of a spin polarization in time stemming from the interaction of nuclear spins 
with a basic dissipative environment, for example, with lattice vibrations. Interestingly, for- 
mulae for the relaxation rate '~fi{uj) (B19) and for a spectral density Ki{uj) of fluctuations of 
an individual nuclear spin in the time- varying magnetic field of the qubit have been derived 
in the process. These formulae are also useful for understanding the electron spin resonance 
in the flipping magnetic field. For the realistic parameters of the phase qubit based on 
the YBCO crystal [^|6|,[T2| we have estimated a decoherence time that meets the necessary 



condition for the successful operation of the qubit. 
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APPENDIX A: 

In the general case, averaging the expression (6) over an ensemble of all spin config- 
urations, namely, finding the functions (sinQt/Q) and ((1 — cosQt)/Q'^) , can be done by 
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means of Fourier transforms, when the dependence on Pq is transferred from the frequency 
Q2 _ ^2 _|_ ^p2 ^Yie exponent e*^^° , for example, 

, sin Qt , 



with 



-) = {J rfee^«^°F(0) (Al) 



1 /-co sin2t./(A/2)2 + 
F(e) = — / dxcosi^x) ' ' (A2) 



After calculating the integral (A2), we get for F(^) 



^=^(0 = (^V^^^-e) e{2t - lei) (A3) 

where 9{2t—\^\) is the Heaviside step function, and Jq is the Bessel function. In the formula 
(Al) we average the function exp(i^Po) or cos(^Po) over the spin configurations: 

r^)^\[\(^.iP„)J„{^4^). (A4) 
The characteristic functional of the spin ensemble can be represented as 

(e^^^°) = (exp(^eEY(^0.))- 
(exp ^— -(si)^ • ...exp %——[s2n)z \ ) = 



2 ^ ^'^) V 2 

cos ^-^ j . ... cos ^-^ j = n cos J . (A5) 

Here 2N is the total number of the nuclear spins in the qubit's volume. As a result, for the 
averaged function (A4) we obtain 



(A6) 



For the sake of simplicity we restrict ourselves to the case of the homogeneous magnetic 
field created by the flux: B(ri) — ... — B{r2N), and Ai = ... = A2Ar = Aq. It should be 
mentioned that in the real case the dispersion in Aj can be of order of its average value 
(Aj). However, as it will follow from Eq.(lO) derived later in the Section HI, taking into 
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account this dispersion will change an estimation for the dephasing rate Ao; by factor of 
order four that is negligible small compared to possible variation in the number of nuclear 
spins 2N 10^ — 10®. With the above-mentioned approximation we get 

S (¥) = (¥) = p„ ^^'^ - ^^^^ 

with being the binomial coefficients. Calculating of integrals gives us the simple formula 



An expression for ((1 — cosQt)/Q'^) can be derived by the same approach. These formulas 
should be substituted into Eq.(6) in order to find a time evolution of the qubit's variable 
(rAt)). 



APPENDIX B: 

The calculation of the correlation function of i-th spin {{l/2)[{ai)z{t), {ai)z{ti)]+) in the 
presence of flipping magnetic fleld and a heat bath is of interest by itself. This problem could 
be of importance for measuring of the qubit's magnetic field by means of the electronic spin 
resonance (ESR) in the free radical probe coated on the surface of the superconducting 
qubit. The flipping of the qubit's magnetic fleld will be reflected in the characteristics of 
ESR signal. That is why this method can be used for detecting of coherent and incoherent 
oscillations in the qubits. 

To flnd the spectral function oi i — th nuclear spin we start with the Hamiltonian 

Hi = -^Tx{(Ti)z - {(Ti)xQx - {(Ti)yQy - {(Ti)zQz- (Bl) 

The spin operators ((Ji)^, {'^ijy-i (ci)^ obey the following equations 

dx — ^iTxTy — 2(TzQy + 2ayQz, 

dy = -AiT^a^ + 2a^Q^ - 2o-^Q^, 

dz — —2ayQx + 2axQy. (B2) 
18 



Here we omit the index "i" in spin operators. The heat bath operators {Qi)a are supposed to 
be independent for different spins {i = 1, 2, .., 2N] a = x,y, z), so that the response functions 
and the correlators of the free heat bath variables have the form 

imY^Kt), iQ,)fit')U)eit - 1') = ^it,t%A^, (B3) 

((l/2)[(g.)W(t), {Q,)f{t')U) = M{t,m,6^,, (B4) 

M{t,ti) = M{t,ti)e{t-ti),ip{t,t') = (p{t-t'),M{t,t') = M(t-t')- As mentioned above, the 
free heat bath, e.g., phonons without any interaction with spins, is in thermal equilibrium 
at temperature T. In this case the heat bath is characterized by the susceptibility 

as well as the the spectral density of fluctuations 

S{uj) = J dre'^'^MiT) 

which are related by the fluctuation-dissipation theorem 

S{uj) = nx"{uj) coth{nuj/2T). (B5) 

The imaginary part of the susceptibility, x"{^)j describes dissipative properties of the equi- 
librium heat bath. 

Following the method developed in |[16|- |T^] we can derive from Eqs.(B2) the Langevin- 
like equations for the spin operators, especially, for z-projection cr^ we are interested in. Due 
to Gaussian statistics of unperturbed heat bath variables Q^^\t){a = x,y,z) the response 
of the heat bath to the action of the spin is linear in the spin matrices cTq : 

Qa{t) = gi°n^) + / dtiifit - ti)a,(ti). (B6) 

After substituting these operators into Eqs.(B2) we have to eliminate the unperturbed heat 
bath variables Q^^^ (t) to obtain the equations only for spin variables. To do that the prod- 
ucts Q'^^\t)af3{t){P = x,y,z) involved in the Heisenberg equations (B2) is conveniently 
represented as 
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Q''a\t)Mt) = {QrW,^/3W} + / dhM{t,h)zMt),a4h)]_. (B7) 
This formula can be considered as a definition of the brackets {..} having zero mean value: 



because of the quantum Furutsu-Novikov theorem []TB[. Introducing the fluctuation forces, 
for example, 

=2{gf (t),a.(t)}-2{gf (t),a,(t)} (B8) 

with zero mean values, {^z(t)) = 0, we obtain the Langevin-like equation for z-projection of 
i-th spin: 

2 1 dh{M{t,ti)t[a^{t),ay{h)]^ + ip{t,h){l/2)[a^{t),ayih)] + }. (B9) 

The brackets {..} simplify also the problem of calculating correlations functions of the fluctu- 
ation forces. They merely imply that pairings of operators within the same brackets should 
be omitted in the process. As a result, the explicit expression for the fluctuation force C,z 
(B8) together with the definition (B7) allow us to get the following correlator 

{{l/2)[Ut),Ut')U) = 2M{t,t'){[a,it),a,{t')U + [a,(t), (BIO) 

It should be emphasized that this result is not based on the fluctuation-dissipation theorem 
and will be valid in the case of strong non-equilibrium as well. However, here we resort to the 
approximation of weak coupling between the spin an the heat bath. With this approximation 
we are able to calculate the (anti) commutators in Eqs.(B9),(B10) by means of the equations 
describing the spin evolution in the time-varying qubit's magnetic field: 

(^xit) = (^x{ti) cosAi(t,ti) + cryih) sinAi(t,ti), 
ay{t) = (Jy{ti) cosAi(t,ti) - a^{ti) sinAi(t,ti), 

CTzit) = a,(ti), (Bll) 

20 



where 



Ai{t,ti)^Ai f dhT^{t2). 
Ju 



(B12) 



The operators (Bll) obey the equations (B2) where all heat bath variables Qa are omitted. 
Taking into consideration a free oscillations of the qubit operator Tx{t) with the tunneling 

frequency A we obtain for the function Aj(t, ti) : 



A 



h) = [T^{t) sin A(t - ti) + Ty{t){l - cos A(t - h))] . 



(B13) 



With the properties of Pauh matrices we find also that 



cosAj(t,ti) = cos (2-^ sin 



A(t-ti) 



(B14) 



sin Aj(t, ii) = 2 sin ( 2-^ sin 



Tx{t) cos 



A(t-ti; 



+ Ty{t) sin 



A{t-ti 
2 

"A(t-ti 



X 



(B15) 



As a result, for z-projection of i-th spin operator we obtain the simple stochastic equation: 



(B16) 



Here 



Viit) = A j dtiip{t,ti)smAi{t,ti). (B17) 

The relaxation of the i-th nuclear spin is determined by the real part 7i(a;) of the Fourier 
transform of the function 7j(t, ti) : 



7i(a;) = 4 / drMir) cos(a;r) cos 



2-^ sin 

A V 2 y 



After a short calculation we derive the following formula for the function 7i(a;) : 



(B18) 



7^M = 2E'^-(2^ 



b [uj — — \ + b [uj ^ 



(B19) 
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Here Jm{z) is the Bessel function of m-th order, and S{u) is a spectral density of the heat 
bath fluctuations (B5). It should be emphasized that it is the formula (B19) that gives a 
line width of the spin- 1/2 coupled to a heat bath in the presence of the flipping magnetic 
field created by the qubit. This line width can be measured in experiments on electronic or 
nuclear paramagnetic resonance. 

It follows from Eqs.(B10),(Bll) that the Fourier transform of the correlation function of 
the fiuctuation forces is also determined by the formula hke Eq.(B19): 



(B20) 



As a result the spectral density of fiuctuations oii — th nuclear spin 




is determined by the expression 



KM 



whereas the spectral density (15) has the simple form: 



(B21) 




with the function 7i(a;) given by Eq.(B19). 
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